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Abstract 

fSJ ' We demonstrate behaviour of the momentum space screened Coulomb t-matrix, obtained by a numerical 

^ , solution of the three-dimensional Lippmann-Schwinger equation. Examples are given for different types of 

screening. They prove that it is possible to obtain numerically a reliable three-dimensional screened Coulomb 

00 

' t-matrix, what is important in view of its application in few-body calculations. 
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I. INTRODUCTION 



The few nucleon systems can be nowadays studied using realistic nuclear forces in all their 
complexity. However, addition of the Coulomb force acting between protons poses serious obstacle 
and is still a challenging task. Up to now the long range Coulomb potential could be exactly 
implemented only into calculations of '^He and ^He bound states and of scattering states up to 
three nucleons. 

In the proton-proton (pp) elastic scattering Coulomb force can be included exactly and ob- 
servables can be calculated using e.g. the Vincent-Phatak method For three-nucleon (3N) 
reactions it is possible to include the Coulomb force when performing calculations of the elastic 
proton-deuteron (pd) scattering both in a coordinate as well as in a momentum space represen- 
tation 0,3]. The pd breakup observables can be predicted only in the momentum space using two 
approaches Q]. Both rely on the screened Coulomb potential and renormalization, what permits 
to apply standard methods valid for short-range interactions. 

In the formalism of Refs. 

B and 01 the Coulomb interaction enters Faddeev equations via the 
screened Coulomb t-matrix. This t-matrix appears in a partial wave decomposed form as well as 
in its direct three-dimensional form for which no partial wave decomposition is performed. This 
allows to include the full Coulomb interaction, not restricted to the finite number of the lowest 
partial waves. Since that three-dimensional screened Coulomb t-matrix is a basic component of 
that formulation, therefore one has to put a special attention to its numerical realization. 

Precise numerical calculations of the three dimensional screened Coulomb t-matrix are possible 
for any type of screening. They are numerically quite demanding but in principle use the concepts 
developed in [7^ for the Malfliet-Tjion potential. This potential consists of two terms which have the 
same functional form as an exponentially screened Coulomb potential with the screening parameter 
n = 1 (for details see below). Such particular value of the screening parameter n allows one to 
perform many steps analytically when obtaining momentum space partial wave decomposed as 
well as direct three-dimensional t-matrix elements. The analogous procedure can be adapted to 
the higher values of n, as done in [sl, where observables for pp scattering were calculated. Also 
there on-, half- and off-the-energy-shell elements of the t-matrix are compared with the analytical 
expressions in the screening limit from Refs. jgl. [lO|. 

In this paper we present details of three-dimensional screened Coulomb t-matrix calculations. 
The behaviour of the numerically obtained t-matrix and its dependence on the screening parameters 
for different types of screening will be presented. In Sect|lT]we describe the numerical procedure 
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to get the three dimensional screened Coulomb t-matrix. In Sect JIIII the behaviour of resulting 
t-matrices obtained with different forms of screening will be presented and t-matrices will be 
compared to the screened Coulomb potential. In Sect llVI the convergence of the partial wave 
decomposition of the screened Coulomb t-matrix will be shown by comparing it with the generated 
three-dimensional screened Coulomb t-matrix. Finally, we summarize in SectlVl 



II. NUMERICAL CALCULATION OF THE SCREENED COULOMB T-MATRIX. 



The three-dimensional momentum space matrix elements of the screened Coulomb potential 
< p '\V^\p >= V^{p',p,x = cos{9) = p ■ p') and the 3-dimensional screened_^C^ulomb t-matrix 

2 

elements < p ^\t^{E)\p >= t^{p\p^p • p'\ E) at energy E = are related by 

-j^ poo pi 

t^ip',p,x;E) = —v^{p',p,x,l) + / dp"p"^ / dx 

(1) 

with 

v^(j)',p,x',x) = / d(l)V^^{p',p,x'x + ^/l- x''^\J\ - x'^cos^)) . (2) 
h 

We solve Eq. ([TJ by generating its Neumann series and applying Fade summation P] for three 
types of screening: 

1. the sharp cut-off screening 

y/(r) = G(i?-r)^ (3) 
with the unit step function G. This screening leads to matrix elements 

{p,P,x',x) = 1^^^ #(1 - cos{qR)) (4) 



where q = \J p^ + p'"^ — 2pp'{x'x + Vl — x'^Vl — x^cos(p). 

I 1 

2. the localized screening [11], where the transition from the pure Coulomb potential to zero 
values takes place smoothly in a finite interval [R, 3R] 

V^R^r) = '-{e{R - r) + ie(r - R)e{3R - r)(l + sin(^))) . (5) 
r 2 2R 

The corresponding v^{p' ,p,x' ,x) is 

^ e2 27r^ - 8q^R^ - 7r^{cos{qR) + cos{3qR)) 
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3. the exponential screening, dependent on two parameters: the power n and the screening 
radius R: 

Vfir) = '-e-^iir . (7) 
r 

In this case 

v^ip ,p,x ,x) = d4) dr e ^fl^ = / ^4) , (8) 

where the function Infiio) = /o°° dr sm.{qr)e~^R'' 

For the value of n = 1 the integration over r can be performed analytically resulting in 

v^{p',p,x',x) = , (9) 

^ y + p2 - 2pj/x + -^Y - 4p'V(l - - x) 

For n > 1 a two-dimensional numerical integration is required to obtain v^{p' ,p,x' ,x). Due 
to strong oscillations of the integrand in ([8]) the big number of integration r-points is needed to 
achieve sufficient precision. This significantly increases the computer time needed for the t-matrix 
calculation, which has to be done on a big grid of p,p',x and x' points. Typically, we solve the 
Lippmann-Schwinger equation ([1]) using 120 p-points and 190 x-points for the sharply cut off 
potential and 95 p-points with 130 x-points for other screenings, what requires over 1.5 x 10^ 
calculations of the v^{p',p,x' ,x) function. The integration over (j) in ([8|) can be performed with 
relatively small number of (/)-points and thus the whole numerical difficulty is shifted to calculation 
of Infi{q). In order to speed it up we use the following method: in the first step we prepare the 
In,R{c[) on a grid of 300 q-points in the range of 0-100 fm~^. In order to calculate the integral 
over r we use the Filon's integration formula [12] which is dedicated to integrals of the product of 
the sine (or cosine) with some nonoscillatory function f{x). The upper limit of integration r^ax 
is chosen sufficiently large so that the integrand approaches zero (e~^ ^ = lO^^*'). Since the 
resulting function In,R{Q) undergoes changes of 10 orders of magnitude in a rather small region 
of q, it is very difficult to handle it properly in further interpolations and integrations. A way 
out is to perform interpolations for the ratio I'^ji'{q) = In.R.{(l) / Ii,r{(i) with analytically known 
Ii,R = qij^pi^-i ■ Variation of that ratio is much more restricted, as shown in FigHJ and we 

use its polynomial representation to get Il^ji'iq) at any value of q. For each value of n and R we 
divide the interpolation region into some optimal number of intervals, optimizing their length as 
well as degree of the polynomial. Typically we have 6 intervals while the degree of the polynomial 
varies between 6 to 12. This allows us to describe the oscillating function with a sufficiently 

high precision, what is exemplified in Fig. [T]for n=3 and R=120 fm. In addition to the solid line 



for -^3*120 (9) ) in Fig(T]show the -^3120(9) values obtained from the polynomial formula. The 

agreement is perfect. The examples of the polynomial representation of I^^°{q) for two sets of 
screening parameters are given in the Appendix. The important and useful feature of 
and thus also of the fit parameters, is its independence from the scattering energy. Due to that 
the interpolation can be done once for a given set of n and R values. 

Once the I'^ji'{q) is calculated also /n. /?,(<?) is known and the final integration over (j) leads to 
v^{p',p, x' , x). In our calculations we use typically a set of 500 gaussian points for the (/)-integration 
for all types of screening. 

The properties of the resulting screened t-matrices are presented in the next Section. 

III. THE SCREENED COULOMB T-MATRIX PROPERTIES 

In this section the three-dimensional screened Coulomb t-matrix t{p',p,x = cos{6)) will be 
shown as a function of momenta p,p' at given scattering angle 9. We choose as examples of 
backward, intermediate and forward angles the following values of 9: 134°, 45°, 10° and 5°. 

In Fig. [2] the real and imaginary parts of the exponentially screened Coulomb t-matrix with n=4 
and R=20 fm at E=13 MeV are shown. The real part of t has a high and steep maximum at small 
momenta at 9 = 134°, which evolves to a ridge lying along diagonal p' = p for smaller angles. The 
spiky structure seen for the smallest angle comes only from the graphical representation on the finite 
grid of p and p'-points. The increasing range of the ridge shows that action of the screened Coulomb 
force becomes more and more important at bigger momenta when moving to smaller scattering 
angles. The imaginary part has a minimum at the on-shell point p = p' = po (w 0.396 fm^^ for 
E=13 MeV). Its absolute value is about one order of magnitude smaller than the maximum of 
the real part. The minimum of the imaginary part becomes deeper and narrower with decreasing 
angle. 

The similar behaviour exists also for other values of the screening radius R. This is exemplified 
in Fig. [3] for R = 120 fm. Taking higher values of R leads to a more restricted range of momenta, 
at which the real part of t takes large values. However, the maximum of Re{t) is much higher than 
for R = 20 fm. Also the range of momenta, where Im{t) has the deep minimum is much smaller 
than for R = 20 fm. 

In Figs, m and [5] the t-matrix for the sharp cut-off screening at E=13 MeV is shown for two values 
of cut-off parameter: R=20 fm and R=80 fm, respectively. The general picture is similar to that for 
the exponential screening, however, some oscillations are visible for both real and imaginary parts 
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of the t-matrix. For the real part they are clearly visible at small momenta. For the imaginary 
part at bigger scattering angles they are relatively large, when compared to the absolute minimum 
of Im[t) and decrease with decreasing angles. 

The t-matrices obtained using the localized screening of Eq. ([5]) are shown in Fig. [6] (i? = 9 
fm) and Fig. [7] (i? = 55 fm). These values of the localized screening range parameter R = 9 (55) 
fm correspond roughly to the values of the screening radius R = 20 (120) fm for the exponential 
screening. The resulting t-matrices are very similar to those obtained with the exponential screening 
at the same angles. 

For the exponential screening we investigated the dependence of the screened Coulomb t-matrix 
on the value of the screening parameter n. We found only weak dependence, as exemplified in 
Fig. [8] where t-matrices at 9 = 10° are shown for n = 2 and n = 3. This picture can be further 
supplemented by parts of Fig. [2] (n = 4) and Fig. [J] (sharp cut-off, what corresponds to infinite 
value of n) at the same angle. It is seen that there is only small decreasing of the height of the 
diagonal ridge for the real part of the t-matrix at the smallest momenta. The minimum of the 
imaginary part becomes narrower and deeper with increasing n. 

Finally, we checked, how good is the approximation of the three-dimensional screened Coulomb 
t-matrix by the screened Coulomb potential alone. To that aim we looked at the ratio 

^ ^ tcip'^P^x) - V^{p\p,x) 

In Fig. [9] we show A at different scattering angles, obtained for the exponential screening with 
i?=120 fm and n=4 at -E = 13 MeV. At all angles the real part of A does not exceed 12%. The 
imaginary part of A is below 3% for all p and p' values what emphasizes the smallness of the 
imaginary part and would indicate on validity of approximation by pure screened Coulomb 
potential Vj^. The differences between the values of and Vj^ are biggest for p,p' < po and 
around p = po or p' = po- It means that in those regions of momenta the approximation of the 
screened Coulomb t-matrix by the corresponding screened Coulomb potential is rather poor. 

2 

For the on-shell screened Coulomb t-matrix elements < Pop'\t^{^)\poP > obtained with expo- 
nential screening we show quality of the Born approximation in Figs. [TO] and [TTJ There the screened 
Coulomb potential Vj^{pQ,po, x) together with t-matrix elements are shown as a function of cosine 
of scattering angle 9. The screening parameters are n = 4 and = 20 fm for Fig. [TO] and n = 4 
and R = 120 fm for Fig. [U] 

The real part of the screened Coulomb t-matrix (solid line) and the screened Coulomb potential 
(dotted line) are close to each other, and their ratio does not exceed 4% of Vj^ for the screening 



6 



range i? = 20 fm and increases up to about 10% when the screening radius reaches R = 120 fm. 
In both cases the imaginary part of the t-matrix is much smaller than its real part and becomes 
important only at very small scattering angles. The real part of the t-matrix is in all cases somewhat 
smaller than the corresponding potential matrix element. 

IV. THE THREE-DIMENSIONAL T-MATRIX FROM THE PARTIAL WAVE DECOM- 
POSITION 

In this section we would like to compare the on-shell elements of the three-dimensional screened 
Coulomb t-matrix obtained directly from Eq. ([TJ with its value derived from solutions of the partial 
wave decomposed Lippmann-Schwinger equation < p,l\t\p' >. The latter is given as a sum of 
contributions from different angular momenta I 

Imax o/ I 1 

^ Yl -^Pi{x)<pAt\p',i >, (10) 

1=0 ^ 

where Pi{x) is the Legendre polynomial of the cosine of the scattering angle, x = p ■ p', and Imax 
should be high enough to give convergent result at given energy and screening. The < p, l\t\p' , I > 
partial wave t-matrix element is a solution of the one-dimensional Lippmann-Schwinger equation 
with the screened Coulomb potential driven by its partial wave element: 

2a r n 

< p,l\V^\p' ,1 >= — / dr r ji{pr)e~^R^" ji{p'r) . (11) 
71" Jo 

In Fig(12]we show convergence of the result with respect to Imax for Imax = 3, 5 and 10. The par- 
tial wave generated screened Coulomb t-matrix is compared to the exponentially screened Coulomb 
t-matrix obtained directly from Eq. ([T]) with n = 4 and i? = 20 fm at E=13 MeV. For such rela- 
tively small screening radius R it is sufficient to restrict to partial waves up to Imax = 10 only to 
reproduce the full 3-dimensional t-matrix (thick dotted line). However, the number of partial waves 
needed to reproduce the full three-dimensional t-matrix increases rapidly with increasing screening 
radius. This is exemplified in FigJTSl for R = 120 fm, where partial wave generated results with 
Imax = 10, 20, 30 and 50 are shown together with the full solution. It is clearly seen that even 
taking Imax = 50 is still not enough to describe with sufficient precision the full three-dimensional 
screened Coulomb t-matrix. This shows the importance of the direct three-dimensional solution of 
Eq.([T]), which avoids problems caused by very slow convergence of the partial wave expansion for 
big screening radii. 
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V. SUMMARY 



In this paper we investigated numerically behaviour of the three-dimensional screened Coulomb 
t-matrix using different types of screening. Such t-matrix is solution of the three-dimensional 
Lippmann-Schwinger equation and is an important component of a recently developed novel ap- 
proach to include the pp Coulomb force into the three- nucleon Faddeev calculations [j, 6|. The 
direct application of the three-dimensional screened Coulomb t-matrix decreases substantially num- 
ber of partial waves needed in 3N calculations reducing them to the number of partial waves needed 
to get converged results in case when only nuclear part of the potential is acting. The presented nu- 
merical procedure enables to get precise values of the three-dimensional screened Coulomb t-matrix 
for any form of the screening. 

We have also presented behaviour of the screened Coulomb t-matrix for different types of screen- 
ing and for different values of screening parameters. The resulting t-matrices are similar. Only 
the t-matrix based on the sharply cut off Coulomb potential reveals oscillations, which are absent 
for other types of screening. Those oscillations cause, that numerical requirements on computer 
resources are bigger in that case. The behaviour of the real and imaginary parts of the three- 
dimensional screened Coulomb t-matrix with varying scattering angle is clearly seen in given ex- 
amples. The range of relative momenta, where the real part of t is important rapidly grows with 
decreasing scattering angle for all types of screening. Contrary, the imaginary part of t becomes 
more localized with decreasing scattering angle but the depth of its minimum grows. 

The presented results exemplify, that the three-dimensional screened Coulomb t-matrix can be 
obtained numerically in a reliable way for any type of screening making it a valuable input in 
three-body calculations. 
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VI. APPENDIX 

In the following tables we give values of parameters for the polynomial representations 
of the I^^°{q) = Ylk'^kQ^' '^^ ^^^^ exponential screening and two sets of {n,R) parameters: 
(n = 2, R = 20) in Tab. Uand {n = 4, R = 120) in Tab. Ull The k (first column) gives a power of q 
which goes with the corresponding factor given in next columns for different ranges of q [fm^^]. 
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k 


q<0.4 


0.4<q<0.4825 


0.4825<q<0.98 


0.98<q<1.725 


1.725<q<10.5 


10.5<q<100 





0.5000125677171 


2.400709537574 


3.921200500830 


1.100526033577 


1.015402877247 


1.000334483409 


1 


-0.006479230796918 


-11.71854218076 


-37.98859669416 


-0.2600475901233 


-0.1503793799690E-01 


-0.4682291329186E-04 


2 


167.3502341712 


42.56921618797 


242.8850249646 


0.2987640252234 


0.6468810517976E-02 


0.2854148534117E-05 


3 


-20.31247086770 


-80.71301847419 


-974.8378775605 


-0.1809241210266 


-0.1511703363356E-02 


-0.9406333859035E-07 


4 


-12352.87853580 


78.54728957584 


2694.796660079 


0.5644212648115E-01 


0.2040113688732E-03 


0.1787524617018E-08 


5 


36456.32220409 


-31.11252021404 


-5362.934658952 


-0.7169472025225E-02 


-0.1 5845939471 54E-04 


-0.1954522705295E-10 


6 


-524160.9711593 


- 


7842.047876700 


- 


0.6567867383073E-06 


0.1140554415044E-12 


7 


15708513.64668 


- 


-8464.239235675 


- 


-0.1124054147047E-07 


-0.2748652734470E-15 


8 


-165490836.2729 


- 


6679.063924791 


- 


- 


- 


9 


931536159.6442 


- 


-3752.307134393 


- 


- 


- 


10 


-3219232757.359 




1423.137799525 








11 


7104191720.650 




-326.9144001002 








12 


-9827857786.235 




34.37700645114 








13 


7804143302.566 












14 


-2723712664.926 













TABLE I: Parameters for the polynomial representation of I,-^ji°{q) for exponential screening with n = 2 
and R = 20. 



k 


Q<0.0735 


0.0735<q<0.11 


0.11<q<0.148 


0.148<q<0.4755 


0.4755<q<1.49 


1.49<q<100 





0.4430751038110 


-61.12120265116 


17222.40652970 


1.014501500770 


1.002164118242 


1.000039812743 


1 


0.1624327164522 


3122.555537517 


-1665157.950661 


-0.1499375990712 


-0.8401588857950E-02 


-0.1260392857265E-04 


2 


5657.908456700 


-61815.86410617 


72913086.23221 


0.6806385287457 


0.1422552507692E-01 


0.1308082855220E-05 


3 


37890.24852747 


603243.9054240 


-1912425443.395 


-1.605992271344 


-0.1246326266740E-01 


-0.6209139096908E-07 


4 


-14459376.10638 


-2904646.753554 


33474324944.51 


1.925075383945 


0.5526732459266E-02 


0. 1534201 149986E-08 


5 


613454088.3275 


5524183.167489 


-412042007272.1 


-0.9274592364796 


-0.9817848874961E-03 


-0.2041445354694E-10 


6 


-34028042154.54 




3658148517742. 






0.1386086594813E-12 


7 


1709808009779. 




-0.2360571796460E+14 






-0.3765737840762E-15 


8 


-0.5319386618092E-I-14 




0.1098886587674E-fl5 








9 


0.1051045675524E-(-16 




-0.3598714293204E-fl5 








10 


-0.1399718714098E-I-17 




0.7868023497509E-I-15 








11 


0.1294364258655E-(-18 




-0.1030692803221E+16 








12 


-0.8173922266757E-I-18 




0.6113583130648E-I-15 








13 


0.3207993078724E-(-19 












14 


-0.5891645069211E-I-19 













TABLE IL The same as in Tab U but for n = 4 and R = 120. 
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q [fm"'] 

FIG. 1: (color online) The for different combinations of screening parameters and the quality of 

polynomial fit for I^'f^oil)- The solid (black) line represents I^'^^ioiQ) calculated by direct integration over r 
in ([U as a function of momenta q. The black x-es show the values of I^'iioi'l) obtained from its polynomial 
representation. Other curves: dash-dotted (green), double-dashed-dotted (red), dotted (blue) and dashed 
(magenta) represents /a^f 20 > ^4,^20 > ^3,60° and respectively. For 0.15 < g < 100 fm the is 

practically equal 1. 
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Re t(p',p,x) [fm'] 



= 134° 




Im t(p',p,x) [fm^ 
O.Oe+00 
-5.0e-05 
-1.0e-04 
-1.5e-04 
-2.0e-04 
-2.5e-04 
-3.0e-04 
-3.5e-04 



■^.8 

''0.6 
'S-'* p' [fm"'] 



p[fm 



Re t(p',p,x) [fm' 




Im t(p',p,x) [fm ] 
O.Oe+00 





FIG. 2: The real (left) and imaginary (right) parts of the screened Coulomb t-matrix t{p,p',x = cos{6)) at 
E=13 MeV for different scattering angles 6 = 134°(l-st row), 6 = 45°(2-nd row), 9 = 10°(3-rd row) and 
6 = 5° (4-th row). The exponential screening with R=20 fm and n=4 has been applied. 
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Re t(p',p,x) [fm'] 




Re t(p',p,x) [fm'] 




Im t(p',p,x) [fmT 
O.Oe+00 r 




0.8 

'0.6 , 
'0.4 P'[fm 1 



Re t(p',p,x) [fmT 




0.6 

'0.4 p' [fm"'] 



Im t(p',p,x) [fm' 




Im t(p',p,x) [f 





FIG. 3: The same as in FigElbut for R=120 fm. 
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Re t(p',p,x) [fm'l 





Im t(p',p,x) rfm'^l 



Re t(p',p,x) [fm'l 




0.6 

'0.4 p' [fm"'] 




P' [fm-'l 



2^0 





Im t(p',p,x) [f 




0.8 

'0.6 , 
'0.4 P'[fm 1 



p[fm-'] 



FIG. 4: The same as in Figl2]but for the sharp cut-off screening of Eq.© with R=20 fm at E=13 MeV. 
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Re t(p',p,x) [fm' 




Im t(p',p,x) [fm' 




Re t(p',p,x) [fm'] 




Re t(p',p,x) [fmT 




or 



Im t(p',p,x) [fm' 




lmt(p',| 




Re t(p',p,x) [fm' 




0.8 

'0.6 , 
'OA P'[fm"'l 



Im t(p',p,x) [f 




FIG. 5: The same as in Figl2]but for the sharp cut-off screening of Eq.© with R=80 fm at E=13 MeV. 
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FIG. 6: The same as in Figl2]but for the locahzed screemng of Eq.® with R=9 fm. 
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Re t(p',p,x) [fm'] 




Im t(p',p,x) [fm'l 










FIG. 7: The same as in Figl2]but for the localized screening of Eq.([5]) with R=55 fm. 
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FIG. 8: The real and imaginary parts of the exponentially screened Coulomb t-matrix at E=13 MeV for the 
scattering angle ^=10°, R=20 fm and for two values of the screening parameter n = 2 (up) and n = 3 (down). 
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FIG. 9: The real (left) and imaginary (right) parts of the ratio A = {t—v)/v (see text) for different scattering 
angles 6: 134° (1-st row), 45° (2-nd row), 10° (3-rd row) and 5° (4-th row) for the exponential screening with 
R=20 fm and n=4 at E=13 MeV. 
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FIG. 10: The real (solid line) and imaginary (dashed line) parts of the on-shell exponentially screened 
Coulomb t-matrix t^{po,po,x) at E=13 MeV as a function of the cosine of scattering angle 9 {x = cos{9)) 
(upper part). The dotted line is the corresponding screened Coulomb potential Vj^{po,po, x). In the lower 
part of figure the ratio Re{t^) /Vj^ is shown. The parameters of the screening are R=20 fm and n=4. 
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FIG. 11: The same as in FigHOlbut for R=120 fm. 
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FIG. 12: The convergence in the partial waves of the exponentially screened Coulomb t-matrix at E=13 
MeV. The real and imaginary parts of the on-shell elements of the three-dimensional screened Coulomb 
t-matrix are given by thick dotted line. The other lines represent partial wave generated results obtained 
by summation of partial waves up to Imax = 3 (dotted), Imax = 5 (dashed) and Imax = 10 (solid). The 
screening parameters are n = 4 and i? = 20 fm. 
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FIG. 13: The same as in Fig. [T^lbut for R=120 fm. The partial waves generated results are obtained by 
summing up to Imax — 10 (dash-dotted), Imax — 20 (dotted), Imax = 30 (dashed) and Imax = 50 (solid). 
The thick dotted line represents as before the three-dimensional screened Coulomb t-matrix. 
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